Topological Analysis of Linear Polymer Melts 
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We introduce an algorithm for the reduction of computer generated atomistic polymer samples to 
networks of primitive paths. By examining network ensembles of Polyethylene and cis-1,4 Polybu- 
tadiene melts, we quantify the underlying topologies through the radial distribution function of 
entanglements and the distribution of the number of monomers between entanglements. A suitable 
scaling of acquired data leads to a unifying microscopic topological description of both melts. 



Polymer chains cannot cross each other. A successful 
conceptual framework embodying this principle at the 
molecular level is offered by the tube model 0,01 ■ The 
tube model postulates that the mutual uncrossability of 
polymer chains generates topological constraints (TCs), 
referred to as entanglements, which effectively restrict 
individual chain conformations in a curvilinear tube-like 
region enclosing each chain. Lateral chain motion is con- 
fined to the length scale of the tube diameter. Large-scale 
motion is promoted via reptation 0|, an effective one- 
dimensional diffusion of a chain along its tube axis. The 
latter provides a coarse-grained representation that char- 
acterizes the chain topology and is called the Primitive 
Path (PP) . Edwards [3| regarded the PP as the shortest 
path constructed by keeping the chain ends fixed while 
continuously tightening (shrinking) the chain contour, so 
that the resulting path has the same topology relative 
to other chains as the chain itself 0, 1^ . Applying this 
construction for all chains 0, |3| leads to a coarse-grained 
picture of a polymer melt that uncovers its large-scale 
topological substructure. The latter is conceived as a net- 
work of entangled PPs underlying the melt structure 0. 
The tube diameter d is assumed uniform and corresponds 
to the mesh length of this network 

Entanglements dominate the rheological and dynam- 
ical properties of large-molecular weight polymer melts 
0, as well as the ultimate mechanical properties in the 
glassy state of these systems 0. Their topological na- 
ture makes their direct experimental study difficult and 
their microscopic definition elusive In this letter, 
we present microscopically determined distributions that 
describe in a statistical manner the topological state of 
fiexible polymer melts. We utilize an algorithm, referred 
to as CReTA (Contour Reduction Topological Analysis), 
which is capable of reducing the atomistic configuration 
of a computational polymer sample to a network of corre- 
sponding PPs. Nodal points of the reduced network cor- 
respond to TCs experienced by individual chains, while 
network connectivity defines the underlying melt topol- 
ogy. Topological measures extracted from networks of 
entangled Polyethylene (PE) and cis-1,4 Polybutadiene 
(PB) melts are compared against corresponding exper- 
imental data. Quantitative agreement establishes these 



networks as meaningful, well defined, structural repre- 
sentations of the underlying melt topologies. 

Mesoscopic views where a polymer melt is represented 
as an entanglement network (rubber analogy) have long 
been used as conceptual abstractions for the development 
of phenomenological models (transient network theory) , 
and analytical theories (tube model) in polymer physics. 
These views are guided from the rubber-like response 
of melts in viscoelastic experiments (rubbery plateau in 
the relaxation modulus), and the Random Walk (RW) 
statistics underlying the structure of flexible polymers. 
RW-like chain conflgurations, above some length scale 
show a self-similarity, to which universalities and scal- 
ing laws detactable in physical properties can be traced. 
Here, by examining networks generated from well packed 
RW-like chain conflgurations of chemically different poly- 
mer melts, we characterize the melt topology through 
distributions which, when scaled accordingly, unveil a 
unifying topological description. The reported statistical 
properties of entanglement networks constitute a missing 
link for the construction of mesoscopic simulation mod- 
els of polymer melts Q and glasses suitable for the 
prediction of rheological and large-deformation mechan- 
ical properties. One can use PP statistics in reverse- 
engineering fashion to construct entanglement networks 
obeying the mentioned distributions, for any flexible 
polymer. Microscopically determined PPs offer a promis- 
ing basis [3, llll lul for investigating the general picture 
of entanglements invoked by the tube model. 

Following Edwards's perspective and as proposed in 
[1,01, CReTA provides a solution to the following geo- 
metric problem ^llj: 'Given a set of uncrossable curved 
lines in space, reduce continuously their contour lengths 
keeping their ends flxed, until they become sets of recti- 
linear segments (entanglement strands) coming together 
at nodal points (entanglements)'. For each polymer we 
examine an ensemble of statistically independent atom- 
istic samples. The ensembles have been thermodynam- 
ically and topologically equilibrated at all length scales 
by previous Monte Carlo simulations [l^ ll^ employing 
chain-connectivity altering moves which alter the system 
topology. Through CReTA, atomistic conflgurations are 
transformed to networks, which are then analyzed to ex- 



tract topological properties. 

CReTA implements random aligning string moves and 
hardcore interactions. CH^; monomers are treated as 
united-atom hard spheres of diameter a. Chains are rep- 
resented as series of fused spheres, since in atomistic poly- 
mer models the skeletal bond distance (~ 1.5 A) is lower 
than a ('^ 3.5A). Periodic boundary conditions apply in 
all cubic sample directions. A string is defined as a set of 
771 > 1 consecutive chain atoms. On each move, a string 
is randomly chosen and the string atoms are displaced 
to corresponding equidistant points on the straight fine 
segment joining the atoms on either side of the string. 
Chain ends are fixed throughout the process. To avoid 
chain crossing and preserve system topology, moves lead- 
ing to overlaps between any string atom with any other 
atom belonging to a different chain are rejected. This is 
sufficient for preventing crossing when chains consist of 
fused spheres. 777 is randomly chosen from an interval of 
Am ~ 10 width. Accepted moves result in simultane- 
ous reduction of chain contours and progressive shrink- 
age of unentangled loops to straight strands composed of 
fused spheres. Rejected moves stem from either mutually 
blocked chain parts, which will lead to entanglements, 
or from contacts between unentangled loops, which will 
eventually disappear. To accelerate chain tightening we 
increase string sizes by gradually increasing Atti. When 
chain contour lengths are no longer diminishing, chain 
thickness is reduced by decreasing a and the process 
starts anew. This aids in tightening meshed unentan- 
gled loops, which, although temporarily blocked, do not 
represent true TCs. Upon decreasing cr, to preserve the 
fused sphere sequence of chains we place an auxiliary 
atom between successive skeletal atoms that have lost 
contact. Auxiliary atoms are excluded from latter anal- 
ysis. The whole procedure terminates when a predefined 
tT/ ~ 0.5 A is reached. At this point an underlying net- 
work structure of interpenetrating 'zigzag' shaped PPs 
has been revealed Further a reduction to attain 

the infinitely-thin-continuous line limit would be time 
consuming and superfiuous. CReTA differs from recent 
similar approaches. In contour reduction is achieved 
through the minimization of an elastic energy which does 
not lead to well defined networks of 'kinky' PPs where 
nodal points can be specified. Kroger's approach is 
also geometric and capable of generating networks but, 
unlike CReTA, cannot track a monomer in the system 
both in the melt and network configurations. The latter 
feature offers great potential for studying hidden topo- 
logical correlations in space, time, and along the chains, 
in close connection with the melt structure and dynam- 
ics. The P{n) distribution presented later is an example 
of this. 

As shown in Fig. the resulting PPs consist of short 
pairwise blocked chain sections connected by straight 
strands of fused atoms. By examining blocked regions we 
resolve PP 'contacts' to pairs of neighboring monomers. 
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Figure 1: (color online). PP network view of a PE melt. 

These specify pairwise associated TCs in the monomer se- 
quence of their parent chains. In addition, they represent 
effective spatial localization points of the TCs each chain 
is subjected to. By convention, we refer to these points 
as entanglements. Each PP is characterized by its ends 
and a set of entanglements with specific spatial locations, 
monomer sequence indices, and pairwise associations to 
other entanglements. The whole structure is reduced to a 
network, with nodal points the entanglements, and edges 
the joining entanglement strands (ESs). Topological anal- 
ysis reduces then to an examination of network proper- 
ties. Entanglement spacing (network mesh) on the PP 
contour is deduced both in monomer and length units by 
calculating TVes and dEs, the ensemble average number 
of monomers and end-to-end length, respectively, of an 
ES. Experimentally, entanglement spacing is mainly in- 
ferred either from the ES monomer length TVe, measured 
from the plateau regime of the dynamic shear modulus 
[1^ . or from the tube diameter d (entanglement distance) 
measured as an intermediate dynamic length in neutron 
spin echo studies of the dynamic structure factor 
These are interpreted on the basis of the tube model, 
where PP conformations are considered as RWs and an 
ES is identified with the PP Kuhn segment 0- Here, the 
relevant quantities are the PP Kuhn length d = /L^ 
and the number of monomers in a PP Kuhn segment 
Ne = NR^/L"^. N denotes the chain average number 
of monomers. R^ and L are the PP squared end-to-end 
distance and contour length, respectively, calculated here 
as ensemble averages. By this 'Kuhn mapping' approach, 
entanglement spacing emerges as the mesh spacing of an 
assumed uniform PP network of interpenetrating random 
walks. 

Table n summarizes our results for A^es, Ne, c?es, d, 
the packing length p and the density of the atomistic 
samples, along with corresponding experimental values. 
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Table I: Entanglement spacing in monomer (A''es, Ne), and 
length (dES, d) units of PP networks, and packing length p 
of the initial atomistic samples (equilibrated at T = 450K 
(PE), T = 413K (PB)). IVes, dss refer to the natural net- 
work mesh. Ne, d refer to the PP Kuhn mapping approach. 
Numbers in parentheses are experimental values from [l^ . 
unless otherwise indicated. The calculated polymer densi- 
ties are 0.776(0.766)", 0.778(0.766)'', 0.867(0.826) g/cm^; the 
polydispersity indices are 1.083, 1.000, 1.053; and the num- 
ber of chains for each system studied are 16, 8, 24 (top to 
bottom). 
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p, an effective chain thickness that controls coil pack- 
ing, is defined as l/{pchR^), where Pch is the number 
density of chains. Volumetric, structural and confor- 
mational predictions for both polymers fl^ Q| are in 
excellent agreement with available experimental data. 
From the additional p data presented here, we conclude 
that the interplay between large-scale chain conformation 
and monomer packing is nicely captured by the atom- 
istic ensembles. Turning to networks, PE samples of 
N = 500, 1000 display practically the same mesh val- 
ues, showing that well entangled samples possess a quan- 
titatively similar underlying topology, independently of 
chain length. The small differences are attributed to 
different packing lengths. A^es, ^es, the network nat- 
ural mesh spacing quantities, are much smaller than the 
corresponding PP Kuhn segment quantities N^, d, with 
N-ES OANe for both materials. This happens because 
of directional correlations between ESs in the same PP 
that decay exponentially with ES separation. That is, PP 
conformations are not RWs, and PP Kuhn segments are 
not free of TCs. At the Kuhn scale, where PPs become 
RWs, Ne, d for both polymers are in good agreement with 
experimental data. If we consider our systems as rub- 
ber networks, however, the small value of Nes implies 
a larger plateau modulus Go than found experimentally. 
Interestingly, in recent Brownian dynamics simula- 
tions of 3D PP networks, the linear viscoelastic response 
of PB solutions is reproduced quantitatively only when 
an Ne of about half the experimental value is utilized, as 
found here for N^s- The discrepancy is attributed 6| to 
the affine as opposed to phantom network model experi- 
mental estimate of Go- In this respect, our results are in 
favor of the phantom model, which for a tetrafunctional 
network would estimate approximately half an N^ for the 
same Go 0. 

In Fig. [2 we present the normalized distribution P(n) 



Figure 2: Normalized distribution of the reduced monomer 
distance between entanglements for the polymers of Table Q 



of the reduced number of monomers N^s in an ES, 
n = N-ES I Nes ■ The data from both materials superim- 
pose on each other, falling on a master curve of the form 
P{n) = ^(e~*" - e-^"), with fitted values b = 1.30, 
c = 3.78. In view of the large difference in iVEs,P between 
PE and PB, the collapse of data suggests a universal char- 
acter of P{n) for linear polymers. The distribution is 
very broad with an exponential tail and in contrast with 
the tube model implies a non-uniform network mesh and 
large fluctuations in the number of monomers of an ES. 
The origin of such fluctuations Q| and their effect on The- 
ological properties have recently been examined (tRIT^. 
Schieber considered a Gaussian chain with a con- 
stant number of Kuhn segments in contact with a 'bath' 
of entanglements. By flxing the corresponding chemical 
potential, he came up with an exponential distribution 
for the number of Kuhn segments in an ES. Stochasti- 
cally, this distribution can be generated by marching 
from one chain end to the other and placing entangle- 
ments on Kuhn segments according to a Poisson process 
of 'entanglement events'. In our case the exponential 
tail of P{n) is the result of a similar stochastic process 
evolving on the monomer sequence space of a chain (see 
below). However, in a Poisson process successive events 
can come arbitrarily close, while it is meaningful to as- 
sume that an entanglement requires a certain number of 
successive monomers to develop. This creates an ejfec- 
tive repulsion between successive entanglements in the 
monomer sequence space of a chain, which will become 
prominent close to and below the number of monomers 
of the Kuhn segment of the atomistic chain. This is ev- 
idenced in Fig. Elby the downturn and vanishing trend 
of P{n) data as n — > 0. 
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Figure 3: Total and partial RDFs of PE-1000 (solid) and 
PB-1000 (dashed). Partial RDFs depict correlations between 
different pair types of binary entanglements. RDF-0 curves 
correspond to entanglements that are 'spatial neighbors', i.e., 
all four chains passing through them are different. RDF-1, 
RDF-2 curves correspond to entanglements at any distance 
that are 'topological neighbors', i.e., connected by one or two 
common chains, respectively, passing through them. Each 
case is sketched in the figure. 

A stochastic interpretation of our results can be given 
in terms of a renewal process , which is the general- 
ization of a Poisson process. In such a process, succes- 
sive events are still identically, independently distributed, 
though with a distribution that is not exponential, cf. 
the normalized P{n) presented here. Moreover, P{n) is 
the convolution of two exponential distributions, namely 
foe"''" and ce~'^"', and can be interpreted as the re- 
sult of two (uncorrelated) alternating Poisson processes 
with rates b, c, evolving on the monomer sequence space 
of a chain. The first process, with rate c, does not result 
in observable events What it does is to stochasti- 

cally create an unentangled monomer sequence in front 
of an entanglement, corresponding to the mentioned re- 
pulsion. The second process, with rate b, takes over when 
an unentangled sequence has been created and places an 
entanglement in one of the following monomers. 

Finally, in Fig. El we present the total and partial ra- 
dial distribution functions (RDFs) of entanglements ver- 
sus r' = r/dES, the spatial distance r reduced by the 
average network mesh length. Apart from some minor 
differences, the curves for PE, PB superimpose on each 
other. Therefore, the spatial network nodal correlations 
stemming from the underlying topologies of the PE, PB 
melts are similar when scaled accordingly. The pres- 
ence of very small network strands, which connect first 
topological neighbors (see Fig. 01, is responsible for the 



effective spatial attraction evidenced in the total RDF 
for r' < 0.4, and separately in RDF-1, 2. The fact that 
RDF-0 curves are below RDF-1 curves for r' < 0.7 indi- 
cates that in the region around and near an entanglement 
topological neighbors persist and spatial neighbors can- 
not easily penetrate. The latter populate primarily the 
region beyond r' ~ 0.7 around a network node, where 
RDF-0 curves overpass RDF-1 curves. RDF-2 curves fall 
rapidly with r' and correspond to pairs of chains inter- 
twined through two entanglements. Intertwining results 
in closely spaced entanglements and is highly improba- 
ble at large r'. The total RDF implies that in the range 
r' > 0.8 the network structure can be described as an 
ideal dilute gas of entanglements 0, with short range 
correlations occurring at smaller r' . 

In summary, common organization properties of topo- 
logical constraints in polymer melts have been revealed 
through a statistical analysis of primitive paths. We are 
grateful to N. Karayannis, V. Mavrantzas and P. Gestoso 
for providing the atomistic polymer samples. This work 
is funded by the Dutch Polymer Institute. 
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